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Recent analysis of relativistically expanding shells of cosmological 
7-ray bursts standard and not peak luminosity (Lq). Assuming a flat 
' Friedmann cosmology (q a — 1/2, A = 0) and constant rate density 

^\ . (po) of bursting sources, we fit a standard candle energy to a uniformly 

t-H ■ selected log iV-log S in the BATSE 3B catalog correcting for fluence 

j_( ' efficiency and averaging over 48 observed spectral shapes. We find the 

Ohi data consistent with _Bo = 7.3^^ x 10 51 ergs and discuss implications 

■ of this energy for cosmological models of 7-ray bursts. 

i> ; 

INTRODUCTION 

t-H ■ On the basis of strong threshold effects of detectors, Klebesadel, Fenimore, 
^ \ and Laros (1982) concluded that GRB fluence tests were largely inconclusive. 

■ As a result, nearly all subsequent number- brightness tests have used peak flux 
\ (P) rather than fluence (S). However, the standard candle peak luminosity as- 

■^j- ■ sumption that is required by log AT-log P studies is unphysical. If, for instance, 
I bursts originate at cosmological distances and are produced by colliding neu- 

■ tron stars (|ll|) the n one might expect that total energy would be standard 
I and not peak luminosity. Moreover, recent analysis of the time histories in 

^ ■ relativistically expanding shell models has found the required differences in 
Q-i bulk r factor between different GRBs all but eliminates the possibility of a 
q ■ standard candle luminosity in such models (g). 

In this paper, we seek to eliminate the large threshold effects present in 
c/3 log iV-log S studies by correcting the observed number of bursts at a given 
fluence by the trigger efficiency of the detector. In §1 we use the calculated 
trigger efficiency in PVO (||) and the catalogue of PVO events (Q) to test the 
correction algorithm. In §11 we examine the fitting algorithm to the log A^-log 
S curve in BATSE 3b and find a standard candle energy for cosmological 
03 gamma-ray bursts. In §111 we discuss the implications of such an energy and 
the distances implied by the fit. 

I. PVO CONSISTENCY CHECK 

The Pioneer Venus Orbiter (PVO) had a peak flux trigger system (sampled 
on 0.25, 1.0, and 4.0 sec timescales) and was sensitive to bursts down to fluxes 
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of 5 x 10~ 6 erg cm -2 . Despite a substantially lower fluence trigger sensitivity 
range, PVO saw hundreds more bright bursts than BATSE due the relatively 
long on-time and large sky-coverage of PVO. As the bright region of the 
BATSE log iV-log S curve seems to fit a -3/2 power law well, we would expect 
that the entire log TV-log S curve of PVO should show a similar behavior. 

Using the PVO trigger efficiency, e(S), from in t' Zand & Fenimore (1996), 
for each burst i with fluence Si in the PVO catalogue we take the expected 
number of bursts with pre-detection fluence to be -/V ijtr uc = N^ b s /e(Si). We 
then compare the derived log N(> 5)-log S curve with a -3/2 power law 
with arbitrary ^-intercept as seen in figure (lb). Although at lower fluence 
there appears to be a deviation from -3/2, the fit is good: with a Kolmogorov- 
Smirnov (KS) probability of 40% that the corrected distribution comes from a 
-3/2 power law. We derive this KS statistic by finding the maximum distance 
between the corrected and -3/2 distributions (in linear space) down to S = 
10~ 4 ' 5 erg cm~ 2 , the fluence at which e(S) falls to 50%. By comparing the 
distributions down to lower fluences, we find an even better KS fit. This 
is understandable because as we add more bursts to the distribution, small 
deviations at high fluences contribute less to the KS distance parameter. Note 
that although the fit is acceptable, a -3/2 slope is not necessarily required by 
the corrected data. We thus conclude that the trigger efficiency determination 
algorithm from in t' Zand and Fenimore (1996) is sound, at least in PVO. 

II. Deriving the log N-log S Curve 



A. BATSE Trigger Efficiency 



One subtly worth noting is that BATSE trigger efficiencies are model depen- 
dent, ie. they depend on the choice of E and cosmology, since it is necessary 
to know a priori the true underlying distribution of bursts that passes by the 
detector. 
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FIG. 1 The log TV-log S curve 
for PVO. a) The uncorrected 
curve for 293 events in the PVO 
catalogue which shows signif- 
icant departure (KS ~ 0.0) 
from the -3/2 power law (shown 
as solid line) expected from 
BATSE observations, b) The 
corrected curve (KS = 0.40 at 
e{S) = 0.50) using the PVO 
trigger efficiency from in t' Zand 
& Fenimore (1996). Any de- 
viation from a -3/2 power-law 
at low flucnces we attribute to 
an incomplete understanding of 
the trigger efficiency. 
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In fact, the derivation of the PVO e{S) assumed an underlying -3/2 distri- 
bution. Petrosian and Lee (1996) have constructed trigger efficiencies us- 
ing bivariate correlation. While this method does not assume a particular 
cosmology, it does require that GRB brightness and duration are inherently 
uncoupled. Our method does not have this requirement and we make no as- 
sumptions about the bursts other then they are cosmological in origin. The 
BATSE trigger efficiencies could be calculated for any E (go = 1/2, A = 0) 
and two are depicted in figure (2a). Note that the efficiency is nearly unity 
for the several orders of magnitude in fluence. The corrected log A-log S 
curve for BATSE is depicted in figure (2b) for two values of E . Interestingly, 
the two distributions arc nearly identical for most of the fluence range. In 
addition, it is clear that the bend from -3/2 in log A-log S is true. 

B. Standard Candle Energy Fits 

The observed fluence of a source depends strongly on the spectrum, and 
since the observed spectral shape depends on the distance to the object, the 
intrinsic spectrum of a GRB object must be used. In addition, the normal- 
ization and the spectral shape vary over the duration of the burst, adding to 
the uncertainty in analysis. 

Following a similar analysis as in Fenimore and Bloom (1995), we take as 
our baseline spectra averages over the GRB spectra fit by Band ct al. (1993). 
Each such baseline burst has associated with it an observed fluence, Si, and 
an observed spectral shape, (f>i(E). 



FUerce Trigger Effieciency in BATSE 



Corrected Log N-Log S 



Channels 2+3 

256 or 1024 ms Triggei 






1 1 1 1 


: Uncorrected**^. \ 


-3/2 : 


E„ = 1e52 erq 
I ...No Z Effects 




. BATSE 3B 

' I .1 T T 1 1 




-7 -6 -5 -4 
Log fluence (ergs cm" 1 } 



Observed Fluence (ergs cm"') 

FIG. 2: a) The BATSE trigger efficiency for different E and b) the corrected 
fog A-log S curve for 830 BATSE bursts. Along with the uncorrected log TV- 
log S curve, we depict a corrected curve corresponding to an assume standard 
candle energy of Eq of 10 52 ergs (dot-dash line) and a corrected curve where 
the effect of rcdshift on the baseline spectra is removed (dash line) . 

Since each Band et al. (1993) burst spectrum is averaged over the duration 
of the burst, we assume in the following analysis that the spectral shape is 
constant, that is, cj)(E,t s ) ~ N(t s )(f>i(E). The fluences, Si [ergs cm -2 ], are 



available for 48 of the Band et al. (1993) bursts in BATSE 3B (|l0|). 

The observed spectral shape, (f>i(E), will not necessarily come from a burst 
at z ~ especially if Eq is large. Therefore, for a given Eq, Si, and 4>i(E) 
we first solve for the redshifts, Zi, of the baseline events associated with each 
spectral shape. The standard candle energy, Eq, is given by, 



E = AttR 



2000 

l z / N(t s )dt 

J30 



E 



1 + Zi 



dE 



(1) 



where N(t s ) is the normalization of the spectrum (units of ergs keV -1 ) as 
a function of time at the source. The comoving distance, Ri iZ , is defined in 
eq. [2] of Fenimore & Bloom (1995). The energy range used in calculating 
E Q in cq. (0) is taken as 30-2000 keV, since we later compare Eq to standard 
candle peak luminosity found in the same energy band. 

The observed fluence of the ith baseline burst in the energy range 50-300 
keV is given by, 

Si = / N{t obs )dt obs / Efa —^E dE, (2) 

JO J50 L 1 + z i . 

where N(t b s ) is the observed normalization of the spectrum. 

For a given standard candle energy, Eq, we numerically determine the red- 
shift (1 + Zi) of the ith baseline burst using eqs. (|l|, ||) and letting z r = z%. 
Note that (1 + Zi) J N{t s )dt s = J N(t obs )dt obs . 



TABLE 1. Best Fit Distribution of E = 7.0 x 10 51 ergs 
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a In ergs cm 2 

b Bursts with C m i n /C max > 1 on the 256 or 1024 ms timescale in BATSE 3b. 



Instead of assuming a spectral shape at the source, we use an average 
over baseline spectra to compute the number of expected observed bursts, 
AN exp [Sj to Sj + i] in some fluence range [Sj, Sj+i]: 

, Nband ,R(S j+ i) 

AN exp [S j toS j+1 ] = - Yl / 6[^(r)]-^r 2 dr. (3) 

Nband jRQSj) 1 + z r 

where Nban d = 48 is the number of baseline spectra used and p Q is the rate 
density of bursts per comoving volume. The quantity Si(r) is the predicted 
fluence (using eqs. [Q, ^j) of the ith baseline burst if it was at a distance r. 
This distance corresponds to a redshift 1 + z r . 

We construct 11 fluence bins (in BATSE channels 2+3 corresponding to 
approximately 50-300 keV) of roughly equal number of bursts. We select 
bursts with C m i n /C max > 1 on either the 256 or 1024 ms timescale, then 
find a minimized x 2 between the number of predicted bursts and observed 
by varying E . For 9 degrees of freedom we find an acceptable \ 2 — 14.7 
corresponding to a standard candle E = 7.3^^ q x 10 51 ergs. Table (0) gives 
the bin ranges, number of observed bursts per bin, number of predicted bursts 
for the best fit energy, and their implied redshifts. 

III. CONCLUSIONS 

Our fit of E Q = 7.0±°{ 7 x 10 51 [30-2000 keV] ergs seems a plausible number 
on the basis that GRBs last on the average 10 sec and Lq = 4.6 x 10 50 erg s~ x 
from log A- log P studies (||). However, this Eq implies a rather large efficiency 
of energy conversion to 7-rays (~ 10%) if the bursting mechanism is colliding 
neutron stars (M to tai — 2.8M©). Nevertheless, this result would seem to help 
resolve the "no-host" problem (cf. Fenimore 1993 et al.). Interestingly, that 
the dimmest bursts (S ~ 5 x 10 -8 erg cm -2 ) are required to be at a redshift of 
1 + z ~ 6.4 given this Eq, would seem to rule out several cosmological models 
that require GRB progenitors to be within galaxies (although see Lu et al. 
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1996). This surprisingly high redshift is due to the correct blueshifting of the 
baseline spectra back to the source in eq. (1). If we neglect this factor, we 
obtain a smaller, more tenable redshift of the dimmest bursts (1 + z = 5.2). 

Whatever the conclusion about the models, we note two important results. 
First, the bend in the log TV-log S curve in BATSE is real, not an artifact of 
strong threshold effects. This implies that we are seeing cither a truncated 
spatial distribution of GRBs (as in Galactic models) or an effect due to the 
expansion of the universe. The bend might also be caused by a combination of 
rate density or number density evolution, and a study of their possible effects is 
certainly warranted. Secondly, with the availability of Monte Carlo modeling 
of trigger efficiencies, log iV-log S tests need no longer be inconclusive. 
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